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Abstract 

A model describing the evolution of a liquid crystal substance in the nematic 
phase is investigated in terms of three basic state variables: the absolute tem- 
perature , the velocity field u , and the director field d , representing preferred 
orientation of molecules in a neighborhood of any point of a reference domain. 
The time evolution of the velocity field is governed by the incompressible Navier- 
Stokes system, with a non-isotropic stress tensor depending on the gradients of 
the velocity and of the director field d, where the transport (viscosity) coeffi- 
cients vary with temperature. The dynamics of d is described by means of a 
parabolic equation of Ginzburg-Landau type, with a suitable penalization term 
to relax the constraint \d\ = 1 . The system is supplemented by a heat equation, 
where the heat flux is given by a variant of Fourier's law, depending also on the 
director field d. The proposed model is shown compatible with First and Second 
laws of thermodynamics, and the existence of global-in-time weak solutions for 
the resulting PDE system is established, without any essential restriction on the 
size of the data. 
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1 Introduction 



Liquid crystals are materials with rheological properties between a conventional liquid 
and a solid crystal, where large elongate molecules give rise to a preferred orientation. 
Many different types of liquid crystals phases have been observed in practical experi- 
ments, distinguished by their characteristic optical properties (such as birefringence). 
When viewed under a microscope with a polarized light source, different liquid crystal 
phases will appear to have distinct textures. The contrasting areas in the textures 
correspond to domains where the liquid crystals molecules are oriented in different 
directions. Within a specific domain, however, the molecules are well ordered. 

Theoretical studies of these types of materials are motivated by real-world ap- 
plications. Proper functioning of many practical devices relies on optical properties 
of certain liquid crystalline substances in the presence or absence of an electric field. 
Typically, a liquid crystal layer sits between two polarizers that are crossed. The liquid 
crystal alignment is chosen so that its relaxed phase is twisted. This twisted phase 
reorients light passed through the first polarizer, allowing its transmission through the 
second polarizer (and reflecting back to the observer if a reflector is provided). The 
device thus appears transparent. When an electric field is applied to a liquid crystal 
layer, the long molecular axes tend to align parallel to the electric field thus gradually 
untwisting in the center of the liquid crystal layer. In such a state, the liquid crystals 
molecules do not reorient light, so the light polarized at the first polarizer is absorbed 
by the second polarizer, and the device loses transparency with increasing voltage. In 
this way, the electric field can be used to make a pixel switch between transparent or 
opaque on command. Color LCD systems use the same technique, with color filters 
used to generate red, green, and blue pixels. Similar principles can be used to make 
other liquid crystal based optical devices. 

There have been numerous attempts to formulate continuum theories describ- 
ing the behavior of liquid crystals flows. We refer to the seminal papers [U [10], where 
Leslie and Ericksen provide a mathematical description of various properties of these 
materials, and to Lin and Liu [12] for the first attempt to analyze the model math- 
ematically. We point out that, to the present state of knowledge, three main types 
of liquid crystals are distinguished, termed smectic, nematic and cholesteric. The 
smectic phase forms well-defined layers that can slide one over another in a manner 
very similar to that of a soap. The smectics are ordered along one direction. In the 
smectic A phase, the molecules are oriented along the layer normal, while in the smec- 
tic C phase, they are tilted away from the layer normal. These phases are liquid-like 
within the layers. There are many different smectic phases characterized by different 
types and degrees of positional and orientational order. The nematic phase appears 
to be the most common, where the molecules do not exhibit any positional order, but 
they have long-range orientational order. Thus, the molecules flow and their center 
of mass positions are randomly distributed similarly to a liquid, but they all point in 
the same direction (within each specific domain). Most nematics are uniaxial: they 
possess a preferred axis that is longer, with the other two being equivalent (so they can 
be approximately described as cylinders). Some liquid crystals are biaxial nematics, 
meaning that, in addition to orienting their long axis, they can also orient along a 
secondary axis. Crystals in the cholesteric phase exhibit a twisting of the molecules 
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perpendicular to the director, with the molecular axis parallel to the director. The 
main difference between the nematic and cholesteric phases is that the former is in- 
variant with respect to certain reflections while the latter is not. 

In this paper, we consider the range of temperatures typical for the nematic 
phase. As already pointed out, the nematic liquid crystals are composed of rod- 
like molecules, with the long axes of neighboring molecules aligned. Such a kind of 
anisotropic structure may be described by means of a dimensionless unit vector d, 
called director, that represents the direction of preferred orientation of molecules in a 
neighborhood of any point of a reference domain. In many experiments, the samples of 
nematic liquid crystals consist of slowly moving particles, therefore a relevant approach 
might be to study the behavior of director field d alone in the absence of velocity. 
However, the flow velocity u evidently disturbs the alignment of the molecules and 
also the converse is true: a change in the alignment will produce a perturbation of the 
velocity field u. Hence, both d and u are relevant in the dynamics, and, to a certain 
extent, also the changes of the temperature d (internal energy). We introduce a very 
simple non-isothermal model for nematic liquid crystals in the spirit of the simplified 
version of the Leslie- Ericksen model proposed by Lin and Liu [IT], and subsequently 
studied in jZHSJ. 

In the proposed model, the time evolution of the velocity field u is governed 
by the standard incompressible Navier-Stokes system, with a non-isotropic stress ten- 
sor depending on V x u, V x d, where the transport (viscosity) coefficients vary with 
temperature. The dynamics of d is described by means of a parabolic equation 
of Ginzburg-Landau type, with a suitable penalization term to relax the constraint 
\d\ = 1. The system is supplemented by a heat equation, where the heat flux is given 
by a variant of Fourier's law in which the dependence on the director field d is taken 
into account, see Section [2j Although such a model may seem rather naive from the 
point of view of real-world applications, the present system of equations captures the 
essential mathematical features of the problem, and, last but not least, it is compat- 
ible with all underlying physical principles, in particular with First and Second laws 
of thermodynamics (as we shall see in Section 1272]) . 

Our goal is to show existence of global-in-time weak solutions to the system, 
without any essential restriction imposed on the size of initial data. In order to avoid 
problems caused by the interaction of the fluid with a kinematic boundary, we suppose 
the latter is impermeable and perfectly smooth imposing the complete slip boundary 
conditions on the velocity u . The existence of weak solutions to the standard incom- 
pressible Navier-Stokes system was established in the celebrated paper by Leray [9]. 
One of the major open problems is to clarify whether or not the weak solutions also 
satisfy the corresponding total energy balance, more precisely, if the kinetic energy 
of the system dissipates at the rate given by the viscous stress. Even the so-called 
suitable weak solutions introduced by Caffarelli, Kohn and Nirenberg [3] allow for an 
uncontrolled dissipation of the kinetic energy that may not be captured by any term 
appearing in the classical formulation of the problem. Since the loss of kinetic energy in 
any energetically closed system must be compensated by a source term in the internal 
energy balance, the above mentioned problem causes unsurmountable mathematical 
difficulties whenever the equations for the kinetic and internal (heat) energy are sep- 



3 



arated. To avoid this apparent difficulty, we use the idea proposed in [5] replacing the 
heat equation by the total energy balance. Of course, the price to pay is the explicit 
appearance of the pressure in the total energy balance that must be handled by refined 
arguments. Apart from the fact that the resulting system is mathematically tractable, 
such an approach seems much closer to the physical background of the problem, being 
an exact formulation of the First law of thermodynamics. 

Let us finally mention that, with respect to [5], the main difficulty here consists 
in the proof of sufficiently strong estimates on the director field d in order to pass to the 
limit in the approximate problem. In particular, the celebrated Gagliardo-Nirenberg 
inequality is needed in order to control the strongly nonlinear terms containing V x d 
in both the momentum equation and the internal energy balance (cf. equations (I2.12p 
and (I2TT5D below). 

The organization of the paper is as follows. In Section [21 we derive the model 
following the standard physical principles. The weak formulation, together with the 
main result, are stated in Section [31 In Section HI we derive the necessary a priori 
estimates and establish weak compactness of a family of weak solutions subject to a 
priori bounds. Finally, in Section we introduce a family of approximate problems, 
based on Galerkin-type approximations, and construct a weak solution of the system. 



2 Mathematical model 

Suppose that the fluid occupies a bounded spatial domain Q C M 3 , with a sufficiently 
regular boundary. Let g = g(t, x) and u = u(t, x) denote respectively the mass density 
and the velocity in the Eulerian reference system. Accordingly, the mass conservation 
is expressed by means of continuity equation 

<9 t £ + div(H = 0, (2.1) 

where, in addition, the standard incompressibility constraint 

divu = (2.2) 

is relevant in the context of nematic liquid crystals. 

By virtue of Newton's second law, the balance of momentum reads 

d t (gu) +div(£>u®u) = divT + gf , (2.3) 

where T is the Cauchy stress, and f is a given external force. 

Motivated by Lin and Liu |TT] we consider the stress tensor in the form 

T = S - gA(0) (V x d © V x d) - pi , (2.4) 

where p denotes the pressure, and § is the conventional Newtonian viscous stress 
tensor, 

S(0, V x u) = (V x .u + V^u) . (2.5) 

Note that the transport coefficients \x and A are functions of the absolute temperature 
$ , see also a strongly related model by Blesgen pQ . More specifically, \i is the viscosity 
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coefficient assumed always positive, while A denotes the thermal dilatation coefficient 
that is an increasing function of 

We assume that the driving force governing the dynamics of the director d is 
of "gradient type" ddJ , where the potential J is given by 

J(^g,d) = W(d) + -G(^g). (2.6) 

v 

Here G is a regular function of d and g, and W penalizes the deviation of the length 
\d\ from the value 1. W may be a general function that can be written as a sum 
of a convex (possibly non smooth) part, and a smooth, but possibly non-convex one. 
A typical example is W(d) = (\d\ 2 — l) 2 . Consequently, d satisfies the following 
equation 

d t d + u ■ V x d + d d W(d) = -div (gV x d) . (2.7) 

Q 

Finally, in accordance with the First law of thermodynamics, the internal energy 
balance reads 

dt(ge-mt) + div (^e int u) + div q = T : V x u , (2.8) 

where e- mt denotes the internal energy density and q its flux. Following Ericksen's 
model [1], the flux can be taken in the form 

q = - K {&)V x ti - («|| - K ± )(#)d(d • V^), (2.9) 

where k, k\\ — kj_ are positive functions of the temperature. Finally, we take e int = c v $, 
where c v > is the specific heat at constant volume. 

Scaling the last equation to have c„ = 1 , we arrive at the following system: 

dtg + div (gu) = , (2.10) 

divu = 0, (2.11) 

d t (gu) + div (gu ® u) + Vp = divS - div (g\($) (V x d V x d)) + gf , (2.12) 

dt(g&) + div(^u) + divq = § : V x .u - gA(0)(V*d V x .d) : V,u, (2.13) 

a t d + u • V^d + aVT(d) = -div (^V x d) . (2.14) 

Q 

2.1 Boundary conditions 

Equations fl2.10l - l2~T4l) must be supplemented by a suitable set of boundary conditions. 
In order to avoid the occurrence of boundary layers, we suppose complete slip boundary 
conditions for the velocity 

u ■ n| an = 0, [Tn] x n| an = 0, (2.15) 
together with the no-flux boundary condition for the temperature 

q-n| sn = 0, (2.16) 
and the Neumann boundary condition for the director field 

VJi • n| an = for i = 1, 2, 3. (2.17) 
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The last relation accounts for the fact that there is no contribution to the surface 
force Tn from the director d. This type of boundary conditions not only simplifies 
the analysis but it is also suitable for implementation of a numerical scheme (cf. [13] 
for further comments on this topic). 



2.2 Energy, entropy 

Multiplying momentum equation (12.121) by u and adding the resulting expression to 
(I2.13P we deduce the total energy balance in the form 

d t (g Q|u| 2 + i^ ) +div (g Q|u| 2 + tfJ u) + div (pu) + div q (2.18) 

= div (Su) - div \q\(&) (V x d V x d) u j + gf ■ u. 

Moreover, using the boundary conditions (12.151) . (I2.16P we may integrate (I2.18P over 
Q to obtain 

d t [ (lg\ u \ 2 + e&) = [ gfu, (2.19) 



in particular, the total energy is a constant of motion as soon as f = 0. 

Let us denote by A(i?) a primitive of 1/A(i9). Testing (12.131) by 1 / A(i9) and 
(I2.14p by (div (gS7 x d) — gW'(d)), integrating the sum of the resulting equations over 
Q, and using the boundary conditions (I2.17p . together with the equation of continuity 
(EHUD, we get 

%d + u • V x d) (div (gV x d) - gdW(d)) + d t J^qM*)) + J Q " V ^ 

[ -\div{gV x d) - gdW(d)\ 2 + [ S : V*u - [ g{V x dQ V x d) : V^u, 
Jn Q Jn A(i?) Jn 



(2.20) 



and 

f (d t d + u ■ V x d) (div (gV x d) - gdW(d)) 
Jn 

Thus, finally, we arrive at 



dt I g ( A(0) - - ^(d)") = / - |div( ff V x d) - ^W(d)| 2 (2.21) 



12 

2~ c " v ~7 JnQ 
where the quantity 

1 S = f ? (A(^)-|V a; d| 2 /2-^(d)) 

is the entropy density of the system. Accordingly, the expression under the integral 
signs on the right-hand side of (12.211) represents the entropy production. By virtue 
of the Second law of thermodynamics, the entropy production is non-negative for any 
physically admissible process, in particular, we need A' > 0. 
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3 Main results 



For the sake of simplicity, we restrict ourselves to the case of constant density, say 
q = 1, and f = 0. Accordingly, the problem (12.1 11427111) . supplemented with the 
boundary conditions (12.151) . (I2.16p . and the initial conditions 

u(0, •) = u , d(0, •) = d , 0(0, •) = , (3.1) 

reads as follows 

divu = 0, (3.2) 

d t u + div (u ® u) + Vp = div§ - div (A(0) (V x d © V x d)) , (3.3) 

drf + div (0u) + divq = S : V x .u - \($)(V x d V x d) : V x ii , (3.4) 

d t d + u- V x d+dW(d) = Ad, (3.5) 

coupled with the boundary conditions (I2.15H2~TT1) and the initial conditions (13. ip . 

To begin, we introduce a weak formulation of (l3.2H3.5p and formulate our main 
result on the existence of global-in-time weak solutions, without any restriction im- 
posed on the initial data. 

3.1 Weak formulation 

In the weak formulation, momentum equation (I2.3p . with the incompressibility con- 
straint (12.21) . and the boundary conditions (12.151) . are replaced by a family of integral 
identities 

u(t, •) • V x f = for a.a. t G (0, T) (3.6) 



n 



for any test function if G C°°(fl) , 

/ / (u-dt<p + u®u: V x ip) = [ [j:V x <p- fu -f(0,-) (3.7) 
Jo Jn v ' Jo Jn Jn 

for any if G Cq°([0,T) x Q; M 3 ) , ip ■ ii\q^ = 0, where the Cauchy stress tensor T is 
related to the unknowns through the constitutive equation (12.41) . 

Equation (12.71) holds in the strong sense, thanks to the regularity obtained for 
d. More specifically, we have 

d t d + u • V x d + dW{d) = Ad a.e. in (0, T) x SI, V x d t ■ n\ dQ = , i = 1, 2, 3 . (3.8) 

In the weak formulation the total energy balance (12.181) is replaced by 

So In (G' U|2 + 9tlf + (^' U ' 2 + ^) u ' V ^ + q • V *^) ( 3 - 9 ) 



T n 



Tu-V x <p- / -|uo| 2 + O <p(0,-), 
'o Jn Jn \ z J 

for any if G C^°([0, T) x Q) , where T, q obey (12 .4p . (12.91) . and by the entropy inequality 
<9 t 0+div(0u)+divq> § : V a .u- A(0)(V a; d© V^d) : V x u in V'((0, T) x fi) , (3.10) 
with q as in (12. 9p and S as in (12.51) . A weak solution is a triple (u, d, 0) satisfying 
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3.2 Main existence theorem 

Before formulating the main result of this paper, let us state the list of hypotheses 
imposed on the constitutive functions. We assume that 

WeC 2 (R 3 ), W>0, dW(d) • d > for all |d| > D (3.11) 

for a certain D > . 

In addition, the transport coefficients are continuously different iable functions 
of the absolute temperature satisfying 

< /£ < ti(&) < p, < « < («|| - «±)(0) < « for all i? > (3.12) 

for suitable constants 

Finally, let A G C^QO, +oo)) be such that 

A'(tf) > 0, A'(0) > 0, A(0) = 0, A(0) < A for all d > (3.13) 

for a certain A > 0. 

Our main result reads as follows. 

Theorem 3.1. Let ft C R 3 be a bounded domain of class C 2+u for some v > 0. 
Assume that hypotheses 113. 11\ - \3.13\) are satisfied. Finally, let the initial data be such 
that 

u GL 2 (ft;R 3 ), divu = 0, d G L°° PI W 1 ' 2 (Q; R 3 ), 

(3.14) 

$ G ^(ft), essinfn^o > 0. 

Then problem h3.(M3.l0^i possesses a weak solution (u, d, d) in (0, T) x ft 
belonging to the class 

u G L°°(0, T; L 2 (ft; R 3 )) n L 2 (0, T; W 1,2 (fi)), (3.15) 

d G L°°((0, T) x ft; R 3 ) n L°°(0, T; W 1,2 (ft; R 3 )) n L 2 (0, T; W^ 2 ' 2 (ft; R 3 )), (3.16) 

G L oo (0,T;L 1 (ft))nL p (0,T;W^ 1 ' p (ft)), 1 < p < 5/4, > a.e. in(0,T)xft, (3.17) 
with the pressure p , 

pe L 5/3 ((0,T) x ft). (3.18) 
The rest of the paper is devoted to the proof of Theorem 13.11 

4 A priori bounds 

We establish a number of formal a priori estimates. These will assume a rigorous 
character in the framework of the approximation scheme presented in Section [5] below. 

Combining (12. 19p (multiplied by a positive constant K > 0) with (12.211) we 
obtain the total dissipation balance in the form 

jf (f |u| 2 + (K4 - A(0)) + + W(d)) (r, •) (4.1) 
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JQ 



Ad - dW(d)\ 2 + : V,u - q • 

< J (f |u | 2 + (KA - A(0 O )) + + W(d )\ . 

For K sufficiently large, the terms on the left hand side in (14.1 p turn out to be non- 
negative, and, in accordance with hypothesis (I3.14p . the integral on the right-hand 
side is bounded; hence we deduce the a priori bounds 

u G L°°(0, T; L 2 (VL] M 3 )) n L 2 (0, T; W 1 ' 2 ^] M 3 )) n L 10/3 ((0, T) x fi; M 3 ), (4.2) 

0, \og(&)eL°°(0,T;L l (Q)), (4.3) 

d G L°°(0, T; W 1,2 (fi; M 3 )), (4.4) 

where we have used (13.111 - 13 . 1 3[) . 

The next step is to take the scalar product of equation (I3.8P with d yielding 

<9t|d| 2 + u • V x |d| 2 + 2dW(d) • d = A|d| 2 - 2| V^dl 2 . (4.5) 

By virtue of hypothesis (13. lip , we may apply the standard maximum principle to |d| 2 
to obtain 

d G L°°((0,T) x f2;R 3 ). (4.6) 
Now, going back to (14.11) and making use of (14.61) . we get 

d G L 2 (0, T; W 2 ' 2 (Q; M 3 )), (4.7) 

which, together with Gagliardo-Nirenberg interpolation inequality (cf. [T4"l p. 125]) 

||Vd|| L 4 (n) < ci||Ad||^, 2 (n) ||d||J2 (n) + c 2 ||d|| LO o (n) , (4.8) 

gives rise to 

V x d G L 4 ((0,T) x H). (4.9) 

This estimate turns out to be "crucial" in order to obtain a bound for the pressure 
and, in general, for the proof of existence of solutions. 

Thanks to our choice of the slip boundary conditions (12.151) for the velocity, the 
pressure p can be "computed" directly from (13. 8p as the unique solution of the elliptic 
problem 

Ap = div div (s - A(#) V x d © V x d - u <g> uY 
supplemented with the boundary condition 

d n p = (div (§ - A(tf) V x .d V x d - u <g> u)) • n onffi. 

To be precise, the last two relations have to be interpreted in a "very weak" sense. 
Namely, the pressure p is determined through a family of integral identities 

J pAp = J (s - \(-&)V x d V x d - u u) : V 2 v? (4.10) 



for any test function ip G C°°(Q), V x y? ■ h\qq = 0. Consequently, the bounds estab- 
lished in (14.21) . (14.91) may be used, together with the standard elliptic regularity results, 
to conclude that 

p G L 5 / 3 ((0, T) x H). (4.11) 

Finally, we turn attention to the heat equation (I2.13p . Multiplying (I2.13P by 
H'{{)) (for a generic H G C 2 ([0, +oo)) ) we deduce its "renormalized" form 

d t H{$) + div (if(tf)u) + div (#'(tf)q) (4.12) 
+H"($) (*;(tf)|V,tf| 2 + (K|| - K X ){#)\d • V,tf| 2 ) 

= H'(&) (s - \($)V x d V*d) : V x .u in X?'((0, T) x fi). 

The choice Hiti) = (1 + i9) v , rj < 1, in (14.121) . together with the uniform bounds 
obtained in (Q|) . and yield 

V x (l + d) v G L 2 ((0, T) x fi; M 3 ) for any < v < - . (4.13) 

Now, we apply an interpolation argument already exploited in [2j. Using (14.31) and 
fT4~T3]) and interpolating between $ G L°°(0, T; L 1 ^)) and $ v G L^O, T; L 3 (fi)) , for 
v G (0, 1] , we immediately get 

G L 9 ((0, T) x fi) for any 1 < q < 5/3 . (4.14) 

Further, observing that, for all p G [1,5/4) and v > 0, 

/ |W| p <( [ |w| V- 1 ) 2 f / 

J(0,T)xQ \J(0,T)xQ J \J(0,T)xn 

we conclude from (14.131) and (I4.14p that 

G L p ((0,T) x tt;R 3 ) for any 1 < p < 5/4. (4.15) 

The a priori estimates derived in this section coincide with the regularity class 
(l3.15fl3TTgj) . Moreover, it can be shown that the solution set of (I3.6H3.9I) is weakly sta- 
ble (compact) with respect to these bounds, namely, any sequence of (weak) solutions 
that complies with uniform bounds established above has a subsequence that con- 
verges to some limit. Leaving the proof of weak sequential stability to the interested 
reader, we pass directly to the proof of Theorem 13.11 constructing a suitable family of 
approximate problems whose solutions weakly converges (up to subsequences) to limit 
functions which solve the problem in the weak sense specified in Subsection 13. 1[ 

5 Approximations 

Solutions to the Navier-Stokes system (13. 6p . (13. 7p will be constructed by means of 
the nowadays standard Faedo-Galerkin approximation scheme, sec Tcmam [16J. Let 
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W^(fl; M 3 ) be the Sobolev space of solenoidal functions satisfying the impermeability 
boundary condition, specifically, 

W$ = {v e W 1,2 (fi; M 3 ) | divv = a.e. in Q, v • n\ dn = 0} . 

Since dQ is of class C 2+u , there exists an orthonormal basis {v n }^L 1 of the Hilbert 
space W$ such that v n e C 2+u , see [H Theorem 10.13]. We take M < N and denote 
Xtv = span{v„}^ =1 . Our strategy is to pass to the limit first for N — > oo and then 
for M — > oo . 

The aproximate velocity fields Ujv.m £ C 1 ([0, X^r) solve the Faedo-Galerkin 

system 

— / u A r iA/ -v = / u^jw ® [uat,m]m : VjjV (5.1) 



dt 



J /u(i?jv,m) (v^Uiv.M + V^u^m) : V x .v + ^ A(^ iV ,M)V :E djv,M © V^d/v.M : V x v, 



L 

Jn 



Uat,m(0, •) ■ v = / u • v 



for any v G X N . Here, the symbol [v] A/ denotes the orthogonal projection onto the 
space span{K}£f=i- 

The functions djv.M are determined in terms of un,m as the unique solution of 
the parabolic system 

d t d NtM + [u NiM ]m ■ Va.djv.jif + 9W(div,Af) = Adjv.M, (5-2) 
supplemented with 

V x (d NM )i ■ n|an = 0, i = 1, 2, 3, (5.3) 

djv,M(0, •) = do,M, (5-4) 

where d 0i M is a suitable smooth approximation of d . 

Next, given Un,m > djv,M, the temperature $n,m is evaluated my means of the 
heat equation (cf. Ladyzenskaja et al. [HI Chapter V, Theorem 8.1]) 

d t $N,M + div ($n,m u n,m) + divq^M = Sjv.ju : V^u^m (5.5) 

- (A(T?jV,Af)(V x djV,Af © Vjjdjv.iu)) : V^U^a/, 

qjv.M • n|an = 0, (5.6) 
^,a/(0,-) = ^o,m, (5.7) 
where S n ,m = h{$n,m) (V x Ujv,iif + V*.u/v >M ), and 

QjV.Af = -k{$N,m)^x$N,M ~ («|| - «±)(^AT,M)diV i iw(diV,M ' V^jv.m)- 

Finally, the pressure piv,M is determined as the unique solution of a system of 
integral identities 

J Pn,mA<P = J (^NM-K^NM)^xd N ,M&V x d NM -u NM ®[n NM ] M ^ : V^.y? (5.8) 
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satisfied for any test function <p G C°°(ft) , V x cp ■ u\qq = 0. In particular, we immedi- 
ately deduce the estimate 

\\PN,M\\L2((0,T)xn) < C(M) . 

Now, taking H(&) = (1 + tf) 1 ', with v G (0, 1/2), in (12712]) . we get 

\\9t^N,M\\(C (lO,T];W^-(n))r < Cll^t^Mllil^CTJxn) < C, 

where C is a positive constant independent of A" and M, with r G (3,+oo), v G 
(0,1/2). 

Regularizing the convective terms in (15. ip . (15. 2p is in the spirit of Leray's original 
approach [9] to the Navier-Stokes system. As a result, we recover the internal energy 
equality at the level of the limit A" — > oo. This fact, in turn, enables us to replace 
the internal energy equation (15.51) by the total energy balance before performing the 
limit M — > oo . For fixed M, N , problem (15.11 - 15. 8p can be solved by means of 
a simple fixed point argument, exactly as in [BJ Chapter 3]. Note that all a priori 
bounds derived formally in Section H] apply to our approximate problem. Thus given 
u G C([0,T)-X N ), we find d = d[u] solving fl5j - E3|) , and then d = tf[u,d] and 
the pressure p satisfying (15.51 - [578]) . Plugging these d, °& in (15 .ip we may find a new 
function T[u] defining thus a mapping u i— > T[u] . Given the a priori bounds obtained 
in SectionHl we can easily show that T possesses a fixed point by means of the classical 
Schauder's argument, at least on a possibly short time interval. However, using once 
more the a priori estimates we easily conclude that the approximate solutions can be 
extended on any fixed time interval [0, T] (see [6j Chapter 6] for details). 

5.1 Passage to the limit as N — > oo 

Having established the existence of the approximate solutions Ujv,m, djv,M, $n,m, 
and Pn,m, we let A" — > oo and use the uniform bounds established in Section H] to 
obtain 

utv,a/ -> u M weakly-(*) in L°°(0, T; L 2 (fi; M 3 )) n L 2 (0, T; W 1,2 (fl; M 3 )) , (5.9) 

StUjv.M -> d t u M weakly in L 2 (0,T; (H /1,2 (fi; IR 3 ))*) , (5.10) 

Pn,m -> Pm weakly in L 2 ((0, T) x fi) , (5.11) 

^,a/^m weakly-(*) in L 2 (0, T; H /1,2 (f2)) n L°°(0, T; L l l v (Vt)) , (5.12) 

- weakly-(*) in (0,(0, T; VK 1,r (fi)))* , (5.13) 

d N , M -> d A/ weakly-(*) in L°°(0, T; H/ 1 ' 2 (fi; M 3 )) n L 4 (0, T; H/ 2 ' 4 (fi; M 3 )) , (5.14) 

d t d NM -> 9 t d M weakly in L 4 (0, T; L 4 (fi; M 3 )) , (5.15) 

for any 1/ G (0,1/2), and r > 3. Note that at this stage M remains fixed in the 
convective term Ujv,m ® [uat ; m]m- 

Hence, applying the Aubin-Lions compactness lemma (cf. fl5"| ) . we deduce that 

Ujv,m -> u M strongly in L 2 (0, T\ L 2 (Q; R 3 )) , (5.16) 
$n,m -> # strongly in L p ((0, T) x Q) (5.17) 
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for any p G [1,5/3). Moreover, at this level of approximation, the director field 
is regular, and, in particular, we have 

Va-djv.M -> V x d A/ strongly in L 4 ((0, T) x O) . 

Hence, we can perform the limit passage 

AT, AT 

-> A(^ M )(V :c d M V x .d A/ ) : V,u M in, say, ^((0,7) x fi). 

Thus we may infer that the limit quantities u M , d M , $ M , and pu solve the 
problem 



Ujif(t, •) • V x (f = for a. a. t G (0, T) 
for any test function tp G C°°(Q) , 

/ / (u M ■ + u M [u M ] M : V^y?) = / / T M : - / u • y?(0, ■) 
Jo Jn y ' Jo Jn Jn 

for any </? G C °°([0, T) x fi; M 3 ) , yp • n\ dn = 0, 

d t d M + [u m ]m ■ V s d M + 9W(dAf) = Ad M , a.e. in (0, T) x ft, 
supplemented with 

V x {d M )i ■ n\an = 0, i = 1, 2, 3, 
dA/(0, •) = d ,M, 

and 

<%$ M + div ($ m Um) + div q M > § M : V z u M - A(t? M ) (V^d^ V x d M ) : V x u M 
in the sense of distributions with non-negative test functions, 

■ n\an = 0, 

together with the total energy balance 



^/ n (i|u„r + ^)=0, 



where 

Tm = S M - \(& M ) {V x d M V x d M ) - p M I 

and 

Sm = K&m) (V^ua/ + V* u M ) . 
13 



5.18) 



5.19) 

5.20) 

5.21) 
5.22) 

5.23) 

5.24) 
5.25) 

5.26) 

5.27) 
5.28) 



Moreover, since the convective term u M <g> [u M ] M is regular (and consequently 
d t u M E L 2 (0,T; (W 1 ' 2 ^;^ 3 ))*)), we can take u M as a test function in (15TB to 
recover the kinetic energy balance in the form 



|2 

Il 2 (c) 



\\uM(t)\\ 2 L 2 (Q) + [ f K$m)\V x u m + V^u M | 2 = ||u | 

+ 2 / / A(^ M )(V x d A/ V,d M ) : V x u M . (5.29) 
Jo Jn 

Similarly, taking v = Ujv.m in (15. ip we get 

\\ u N,M(t)\\ 2 L 2^ + / / /i(^M,7v)|V a: U7V iM + V^U^Af | 2 = ||uo|| 2 2( a) 

Jo Jn 

+ 2 / / A(^ A r iA/ )(V x d A r >M V^djv.M) : V^Ujv.M • 

Jo Jn 

Passing to the limit AT — > oo in the last equation, using (15.161) and (13.131) . and com- 
paring the result with (15.291) . we conclude, by means of (13 . 1 2j) . that 

V x u n ,m -> V x u M strongly in L 2 ((0, T) x fi) . (5.30) 

Accordingly, the inequality (I5.23P may be replaced by 

9 i ^A/ + div(^ M UM) + divq A / = S M ■ V x u M - \($ M )(V x d M (D V x d M ) : V x u M ■ (5.31) 

As a matter of fact, (I5T3TD already follows from (l5^3|) . (15^61) . and (15^91) . 

Moreover, taking in (15.191) Um<P (with <£> e ^((0, T) x f2) ) in place of we get 

<9 t Q| u a/| 2 + tfM^ + div ( Q|u A /| 2 + ^Af^ [u m ]m) + div (p M u M ) + divq M (5.32) 

= div (§a/Ua/) - div (\(# M ) (V x d M V x d M ) u M ) in X>'((0, T) x fi) . 
This concludes the passage to the limit for A" — > oo . 



5.2 Passage to the limit as M — > oo 

Our final goal is to let M ->• oo (I5TM15T21)|) . We notice that the limits in ABU . 
(l5TT2ti5TT3|) . (15TT6H5T7I) still hold when letting M -> oo. Moreover, we have 

ftuA/ -> 9 4 u weakly in L 5/3 (0, T; W~ 1,5/3 (fi; M 3 )) , (5.33) 

Pa/ -> V weakly in L 5/3 ((0, T) x fi) , (5.34) 

d M -> d weakly-(*) in L°°((0, T) x fi; R 3 ) n L°°(0, T; W 1 ' 2 (fi; M 3 )) (5.35) 

nL 2 (0,T; W 2fi {Q;R 3 )) , 

d t d M -> <9 t d weakly in L 2 (0, T; L 3/2 (ft; M 3 )) . (5.36) 
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Now, we can easily pass to the limit M — > oo in (j5.18ff5T2"2"|) to recover fl3.6H3.8l) . 
In addition, by virtue of (Oil . (Oil . (Q3jl . we get 



|Um 
2 



2 



Pa/) [um]m j bounded in L 10 / 9 ((0, Tjxfl), 

/ J M>0 

{^ M u A/ } M>0 bounded in L q (0,T; L q ({l)) for any g G [1, 10/9) , 
{§mU M } m>0 bounded in L 5/4 ((0, T) X tt) , 

{\(#M)(V x d M © V x d M )u M } M>0 bounded in L 5 / 4 (0, T; L 5 / 4 (fi)) . 

Consequently, we can pass to the limit in (15 .32 j) to deduce the desired conclusion (13.91) . 

Finally, as convex functionals are weakly lower semicontinuous, we can see that 
(I5.23P gives rise to 

<9 t tf + div(tfu) + divq> S : V«u-A(i?)(V x d0V x d) : V x u in P'((0, T) x fi) . (5.37) 



This completes the proof of Theorem 13.11 
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